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Abstract 

A simple method is proposed to prepare conveniently the effective pure 
state |00...0)(0...00| with any number of qubits in a quantum ensemble. 
The preparation is based on the temporal averaging (Knill, Chuang, and 
Laflamme, Phys. Rev. A 57, 3348 (1998)). The quantum circuit to prepare 
the effective pure state is designed in a unified and systematical form and 
is explicitly decomposed completely into a product of a series of one-qubit 
quantum gates and the two-qubit diagonal quantum gates. The preparation 
could be programmed and implemented conveniently on an NMR quantum 
computer. 

PASC number(s): 03.67.Lx, 76.60.-k 
1. Introduction 

Quantum computers can solve certain problems that can not be solved by 
any classical digital computers [1-5]. Particularly, the prime factorization of a 
large number can be performed in a polynomial time on a quantum computer 
[3]. This has stimulated the experimental work of quantum computation 
in various physical disciplines including quantum optics [6], trapped ions 
[7], nuclear spin system [8], superconducting Josephson junctions [9], and 
so forth. Quantum computation is usually performed in a pure quantum 
state at any time [10]. However, recently it has been shown that quantum 
computation may be implemented in a spin quanutm ensemble at a finite 
temperature by using nuclear magnetic resonance (NMR) techniques [11, 12]. 
After that, there has been a flood of experimental work of ensemble quantum 
computation on few-qubit systems [13-20] . The key of the approach is that a 
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quantum ensemble can be prepared in certain mixed states, i.e., the effective 
pure states [11] or the pseudopure states [12]. An effective pure state is 
a mixed state in a quantum ensemble that behaves for all computational 
purposes as a pure state, that is, the effective pure state is isomorphic to the 
pure state. Therefore, the reversible quantum computation can be performed 
in an effective pure state just like in a pure quantum state. The key step in 
the approach is then the preparation of the effective pure states. There are 
several methods to prepare the effective pure states or the pseudopure states. 
These include logical labeling [11, 21], spatial averaging [12, 22], temporal 
averaging [23], spectral labeling [17], and so forth. The temporal averaging 
method is really a technique much similar to the time averaging, coherent 
accumulation, and phase cycling techniques in NMR experiments [24, 25], 
but used for preparation of the effective pure states. This method has been 
investigated thoroughly by Knill, Chuang, and Lafiamme [23]. The method 
such as the exhaustive averaging method usually involves cyclicly permuting 
the nonground states in 2" — 1 different ways such that the average of the 
prepared states is an effective pure state. Therefore, number of experiments 
to prepare the effective pure state grows exponentially as the qubit numbers 
n, showing that this method may be reasonable to consider implementing 
the preparation for small numbers of qubits. The advantage of the method 
over other preparations [11, 12] may be that it is a simple method to give a 
higher signal-to-noise ratio for the NMR experimental measurement. 

In this paper a simple method is proposed to prepare conveniently the 
effective pure state p^fj = A|00...0)(0...00| in a spin quantum ensemble that 
consists of a large amount of two-state n-spin (1=1/2) quantum systems corre- 
sponding to the pure quantum ground state |00...0) of the quantum system. 
It is also based on the temporal averaging [23]. However, the preparation 
starts generally from the specific reduced density operator of the longitudi- 
nal n-spin order component cr(0) = 2^~^Iizl2z---Inz [26, 27] that can be easily 
obtained from the equilibrium magnetization in the spin quantum ensemble 
(called also the spin system without confusion). The quantum circuit to 
prepare the effective pure state with any number of qubits is designed in a 
unified and systematical form. It is then explicitly decomposed completely 
into a product of a series of onc-qubit quantum gates and the two-qubit 
diagonal gates [26]. Therefore, the preparation could be programmed and 
implemented conveniently on an NMR quantum computer. 

2. Prepeiration of the effective pure states 
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An effective pure state or pseudopure state Pg^j(^) in a quantum en- 
semble corresponding to tfie pure quantum state |^) may be defined as [11, 
12] 

p,ff{^) = {i-x)E+m{^\ (1) 

wfiere tfie normalization factor is omitted without confusion, A is a real 
constant, and E is N {= 2")-dimensional unity operator. Because the unity 
operator is unobservable and keeps unchanged when it is acted on by an 
arbitrary unitary transformation the mixed state A|^)(^| and the effective 
pure state p^fj{^) cannot be distinguished by standard NMR experiments. 
Therefore, the process of NMR ensemble quantum computing can be really 
characterized completely through the reduced effective pure state a^ff — 
A|\I/)(\E'|. In order to make use of the massive quantum parallelism [2] the 
superposition in a quantum system is usually created at the first step in 
quantum computing, while the superposition in a spin quantum ensemble 
can be created conveniently by acting the Walsh-Hadamard transform W [5] 
on the effective pure state corresponding to the pure quantum groundstate 
|00...0) in the quantum system: 

(Te// = A|00...0)(0...00| (2) 
Thus, in the paper a simple method is proposed how to prepare conveniently 
the effective pure state of Eq.(2) in a spin quantum ensemble with any num- 
ber of qubits. 

First of all, the starting density operator p(0) of the spin quantum en- 
semble is subjected to a unitary transformation whose unitary operator 
exp(— iaQ) is constructed with the Hermitian operator Q defined by Qij = 1 
for all indexes i and j : 

1 — cos( (xN^ 

exp(-mQ)p(0)exp(mg) = p(0) -W),Q] + 

sinfaA^), \1 - cos( aN)]^ + sin^ (aN) ^ 
+i^^[p(0), Q]-+^- ^ '-Qpm (3) 

where [p(0), Q]± = p{0)Q±Qp{0) and the following expansion for the unitary 

operator exp{±iaQ) has been introduced in the above transformation, 

/, ^ 1 — expfiiaA^) ^ 
exp(±iQ;g) = E ^ ^-Q (4) 

It follows from the definition of the operator Q, that is, Qij = 1 for all indexes 
i and j, that the product Qp{0)Q can be reduced to the simple form 

Qp{Q)Q = Q/{p(0)} (5) 
where the function /{p(0)} is defined as 
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N 

f{pm = E Pkiio) (6) 

k,l=l 

For convenience, one writes conveniently the density operator as the general 
form 

p{0)^Tr{E}-^E + a{0\ (7) 
where cr(0) is the traceless reduced density operator, Tr{a{0)} = 0. For 
an NMR spin quantum ensemble the high-tcmpcrature approximation holds 
and equation (7) is met. Obviously, it proves easily with the aid of Eq.(7) 
that equation (3) still holds when replacing the reduced density operator p(0) 
with the density operator cr(0). To simphfy further Eq.(3) the operator Q is 
decomposed explicitly as 

Q = N exp{-qFy)Do exp{i^Fy) (8) 

n 

where the operator F„ = E hu, (a* — ^tVi^)-, hu — ko'kLt, cr is the Pauli's 

fc=i 

operator and the diagonal operator Dg is defined as 

Ds = Diag[Q, .., 0, 1, 0, 0], ([D,],, = 1; [Ds]^^ = 0, ? ^ s). 

Then by exploiting Eqs.(5) and (8) one obtains from Eq.(3) that 
[1 - cos{aN)Y + sin^(a7V) ^ r^^g^^ 

= exp(-iiVaL'o)p+(0) exp(iiVaL'o) - /O+(0) 

+ (1 - cos(aiV))[p+(0),Z}o]+ - ^sin(aiV))[p+(0),Do]- (9) 
where p_|_(0) = exp(i|Fy)p(0) exp(— In particular, if the parameter a 
in Eq.(3) is chosen suitably so that cos(q;A'") — —1 and sin(Q;A'") = 0, then 
equation (9) can be further written as with the help of Eq.(7) 

^/{(7(0)}L>o = exp(-iArQ;L»o)(7+(0) eK^iiNaD^) 

-a+(0) + 2[a+(0),Do]+ (10) 
where c7+(0) = exp(i^Fj^)o"(0) ex.Y>{—i^Fy) . Now, the function /{o"(0)} = if 
the density operator (7(0) is taken as an arbitrary operator of the operator 
set G = {Iky, hz, '^hykx, '^hyhy, 21kyliz, 21kzlix, ■■■}, which does not con- 
tain all the x-components and unity operator of the operator set — {E, 
Ikx, '^Ikxhx, ^Ikxhxlmx-i ■■■}■ Then in this case equation (10) can be further 
simplified to the form 

exp(-iiVaZ)o)o-+(0) exp{iNaDQ) - (t+(0) 

= -2[a+(0),Z}o]+ (11) 
When the starting density operator cr(0) is taken as an arbitrary operator of 
set Gx equation (10) becomes an identity although in this case /{cr(0)} ^ 0. 
According to the definition of the diagonal operator Dq and cos(q;A'^) = — 1, 
sin (a A'") = it turns out easily that 
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exp{±iNaDo) = Diag[-1, 1, 1] = -R (12) 
where R is the phase-shift operation defined in the Grover quantum search 
algorithm [5] . By exploiting the definition of diffusion transform D = WRW 
in the Grover quantum search algorithm [5], equation (11) is rewritten as 
(2/iV) [a+(0), NDo\+ = a+{0) - WDWa+{0)WDW (13) 
Equation (13) is the key to prepare conveniently the effective pure state 
p^^^ = A|00...0)(00...0| in the paper. 

It follows from the definition of the operator Dq that the diagonal operator 
Do is really an effective pure state p^j-f = |00...0)(0...00| = Dq. Because Dq 
is a diagonal operator it can be expanded as a sum of the base operators of 
the longitudinal magnetization and spin order (LOMSO) operator subspace 
of the Liouville operator space of the n-spin (1=1/2) system [17, 26, 27] 

NDo E 24,+ E 4IkJiz+ E ShJiJmz + ... (14a) 

fe=l l>k=l m>l>k=l 

= {El + 27i,) <S,{E2 + 2I2,) ... 0(^„ + 2/„,) (146) 
where i^^ is 2 x 2— dimensional unity operator of the kth qubit in the spin sys- 
tem. In order to extract the effective pure state from the term [^^(O), NDq]^ 
on the left-hand side of Eq.(13) the density operator cr+(0) needs to be cho- 
sen suitably. First, a series of the density operators cr+(0) are suitably cho- 
sen, where each density operator a{0) is only taken as an arbitrary operator 
of set G so that /{o"(0)} = can be met, then the term [(x^{0), NDq]^ 
is expanded explicitly by inserting each cr+(0) into the term. By exploit- 
ing these expansions one can construct explicitly the effective pure state 
Peff = A|00...0)(00...0|. The detailed procedure to prepare conveniently and 
explicitly the effective pure state based on equation (13) is illuminated below. 

2.1 The logical labeling effective pure state 

It has been shown that the following state, called the logical labeling 
effective pure state, can be thought of as an effective pure state if one of the 
qubits in the quantum system is used as a logical label [23] 

= + A(|00...0)(0...00| - |11...1)(1...11|). (15) 
It will be shown below that the logical labeling effective pure state of Eq.(15) 
in the n-spin system can be prepared on the basis of Eq.(13). The density 
operator cr+(0) in Eq.(13) is chosen simply as 

(7+(0) = 2^-^hyhy...Iny (16) 

Obviously, the function /{(7(0)} = 0. It is easy to turn out by using Eq.(14) 
that 

[a+(0), iVDo]+ = T'^H-I^I^...!- + 2-\-iYltlt...I+ (17) 
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= 2--\l{I^...I- + Itlt-I^). if n = 4m 

= 2--Hilrl^...I- - Itlt-Ii), if n = 4m + 1 

= -2--\I{I^...I~ + if n = 4m + 2 

= -2-H{I{I,...I^ - /1+/2+.../+), if n = 4m + 3. 
The right-hand side of Eq.(17) is actually an n-qubit maximum entanglement 
state of the system. Therefore, one can prepare conveniently the maximum 
entanglement state through Eq.(13). Now, if n is an even number one makes 
a unitary transformation on the term [a+{0), NDq]^ to convert it into the 
form corresponding to n being an odd number 
exp(-i^F,)(7r72-.../- + Itl^-I^) eMiOF,) = 

^(jrj2-.../--/i+/2+.../+) 

where n6 = 7i/2, that is, exp(±in^) = ±i. Therefore, no matter that n 
is either an odd number or an even number one can prepare further the 
effective pure state by starting from the n-order multiple-quantum coherence 
2"~^i(7f 7^.../^ — I^l2...In)- By making a unitary transformation on the 
multiple-quantum coherence one obtains the logical labeling effective pure 
state of Eq.(15) without the unity operator term 

exp(-if X 2^hJ2....In.)2^-H{I^I^...I- - Itlt-It) 

Xexp(if X 2''Iij2x...Inx) 

= |A^(|00...0)(0...00| - |11...1)(1...11|) (18) 
On the other hand, by using Eq.(13) one can design the quantum circuit 
to prepare conveniently the logical labeling effective pure state of Eq.(18). 
As an example, consider the case of n = 4m + 1, by inserting Eq.(16) into 
Eq.(13) one obtains the logical labelng effective pure state of Eq.(18) 

p,,= (|00...0)(0...00|-|ll...l)(l...ll|) 

= 2«-iJi,J2..../„. -exp(-2f X 2^hJ2....Incc)WDW exp{i^F^) 
x2"-^hJ2z-Inzexp{-i^F,)WDWexp{t^x2^hJ2,...I^,) (19) 
Experimentally, starting from the reduced density operator, i..e., the lon- 
gitudinal n-spin order component a{0) = 2"-~^Iizl2z---Inz that can be pre- 
pared from the thermal equilibrium state in an n-spin quantum ensemble 
one performs two different experiments to prepare the logical labelng ef- 
fective pure state p^^ of Eq.(18) according to Eq.(19), one experiment with 
the identical operation E and another with the unitary operation exp(— x 
2" hxhx - ■ ■Inx)W DW exp(i|Fa;). In comparison with the flip and swap method 
[23] to prepare the logical labeling effective pure state of Eq.(15) the present 
method is still available even when the number of qubits (n) and the polar- 
ization (6) of the quantum system satisfy n5 ~ 1 or when the initial state of 
the system does not have approximate inversion symmetry, but as shown in 



6 



Ref.[23], in these cases the flip and swap method is failure. 

2.2 The effective pure states 

(a) A two-spin system 

The density operator (T+(0) is chosen as (T+(0) = 2Iiyl2z, 2Iizhy, respec- 
tively. One obtains with the help of Eq.(14) 

exp(-if /i,)[2/ij,/2,, NDo]+ exp(2f Ji,) 

= \ND^ -\{Ei- 2Ji,) ®{E2 + 2/2.), (20a) 
fiy.-^{-i\hx)[2hzhy, NDq\+ exp(i|/2a;) 

= \NDo - + 2hz) <S){E2 - 2hz) (206) 
By plus Eqs.(20a) and (20b) and then exploiting Eq.(13) one gets the effective 
pure state of the two-spin system as follows 

p,f^ = 2Do = lE + 2hzhz 

-exp(-2f/i,)l^DW^exp(^f/i,)(2/i,/2,)exp(-^f/i,)iyDiyexp(2f/iJ 
- cxp(-jf /2x.)VrWexp(^f /2.)(2/i,/2,) cxp(-zf /2.)iy Wexp(^f /2.) 
In NMR experiments of ensemble quantum computation the starting reduced 
density operator (t(0) = 2/1^/2^ is first prepared from the equilibrium state in 
the coupled two-spin (1=1/2) system. By starting from the density operator 
(t(0) = 2/12/22, one performs three different experiments to create three dif- 
ferent mixed states, one experiment with identity transformation E and the 
other two with the unitary transformations: exjp{—i\hx)WDW exjp{i^hx) 
and exp(— i|72a:)W^-DVrexp(i|72a;), respectively. Then the created state (t(0) 
= 2Iizl2z minus the other two created states will give the effective pure state 
2Do - \E. 

(b) A three-spin system 

The density operator it+(0) is chosen as o"+(0) = AIiyl2zhzi '^hzhyhzi 
'^hzhzhyi and 4:Iiyl2yl3y, respectively. With the help of Eq.(14) it easily 
turns out for the density operator cr_|_(0) = ilkylizlmz (cyclicly permuting 
k, I, m) that 

exp{-i^lka:)[4:lkylizlmz: N Do]+ ex.p{i^Ikx) 
= |{|(Efc + 2Ikz) 0{Ei + 2Iiz) 0{Em + 2Imz) 
-\{Ek - 2hz) ®{Ei + 2Iiz) ®{Em + 2I^z)} 
= Alkzhzlrnz - Gx:^{-ilhx)W DW ex^{i\Ikx) 

x{41kzlizlmz) exp{-z^Ikx)WDW exp{zpkx) (21a) 
and for the density operator <t+(0) = Alkyliylmy that 

;|exp(-if X SIkxkxImx)WkyIlyImy, A^I^o]+ exp(if X Slkxllxlfux) 
= U-ki^k + 2Ikz) (S>{El + 2Iiz) <^{Em + 2Imz) 
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'^Ikzllzlmz 

-exp(-2f X MkJiJn..)WDWe^^{ilF^) 

x{AhJiJmz) ew{-^lF,)WDW e^il x ShJiJmx) (216) 
where = hx + Iix + Imx and N — 2^. Then exploiting Eqs.(21) the effective 
pure state in the three-spin system can be built up as follows 

p,ff = 22L>o ^\E+ (22 - l)2^h,hJ^, 

X i2^Iul2zhz) exp{-qikx)WDW exp(if 7^,)} 

- exp(-if X 8IiJ2xhx)WDW cxp{~i^F^) 

x{2^hJ2zl3z)cxp{t^F.,)WDWexp{t^ x ShJ^xhx)- 
Experimently, one first prepares the starting reduced density operator o"(0) = 
^^hzhz^Sz from the equilibrium magnetization in a coupled three-spin system. 
Then only five different experiments are performed in order to prepare the ef- 
fective pure state (227:)o-|£;).Ifthe experimental signal-to-noise ratio is high 
enough in the spin system the contribution of component (2^ — 1)2271272273, 
to the effective pure state can be obtained directly by amplifying the ampli- 
tude of the longitudinal three-spin order component 2'^Iizl2zhz- 

(c) A four-spin system 

The density operator o"+(0) is chosen as ct+(0) = SIiyl2zhzhz, 

^hzhyhzhz, ^Iizhzhy^Az, ^Iizhzhzhy, and Slizhyhyhy, ^hyhzhyhy, 

^hyhyhzhy, ^hyhyhyl'iz-i respectively. Then the effective pure state can be 
easily constructed from Eq.(13) with the help of Eq.(14) 
p,ff = 23L>o ^\E + (23 - 1) x 2^hzhzhzhz 

- E {exp(-if7fc,)W^L>iyexp(^|7fc^) 

k=l 

x{Sh-j2zhzhz) cM-^lh-WDW eMi^h-)} 

4 

E {exp(-if X %hxIixImx)W DW exp[-il{Ikx + hx + Imx)] 

m>l>k=l 

X (Shzhzhzhz) exp[if (7fe^ + Ii^ + Imx)]WDW exp(if x Slkxhxlmx)} 
Obviously, to prepare experimentally the effective pure state 237^o — \E in 
the four-spin system one needs to perform nine different experiments directly 
by starting from the density operator a{0) — Sluhzhzhz- 

(d) A five-spin system 

The density operator o'+{0) is chosen as <J+{0) = 2'^Iiyl2zl3zhzhz, 
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2 hzhylszhzhzT'-t 2 hzhzhzhzhy] 2 Iiyl2ylzylizl'6z-----j 2 hzhzhyl'lyhy] 

2^Iiyl2yl3yl4yl5y, lespectively. On the basis of Eqs.(13) and (14) the effective 
pure state can be readily written as 

p,jf = 2^Do = \E + (2^ - 1) X 2''hJ2zhzhzhz 

- E {exp(-2f4,)PyDW^exp(2f4,) 

k=l 

5 

- E {exp(-if X 8IkxIixImx)WDW exp[-if (7fe^ + Ii^ + Imx)] 

m>l>k=l 

x{2'^hzhzhzhzhz) exp[if (7fe^ + Ii^ + 7„^)]W^DW^exp(if x 87fc^7;^7^^)} 
-exp(-if X 327i,72,73,74,75,)W^Z}l^exp(ifF,) 
^{2HiJ2zhzhzhz) exp(-2fF,.)14^7?l^exp(if x m^J^M^^h^ 
This shows that the effective pure state 2'^7^o ~ \E can be prepared by per- 
forming (2^^^ + 1) different experiments in the five-spin system. 

(e) An arbitrary n-spin (n > 6) system 

The density operator o"_|_(0) is chosen as o'_|_(0) = 2"'~^Iiyl2z---Inzi 

on-l T T T T on-l j j j . on-l T T T T T 

^ ^\z^2y^Zz---^nzi ■■■■> ^ ^\z- ■ -^n-lz^ny^ ^ ^\y^2y^Zy^A:Z---^nzf--i 

on— 1 11 1 1 1 . cyn-l 111111 1 

^ J-lz----'-n—Zz-'-n—2y-i-n—ly-'-nyi ^ -'l2/-'2iy-'3j/-'4j/-'5j/-'6^----'n2 v) 

2"'~^lu...ln-5zln-4yln-3yln-2yln-iylny; , respectively. With the aid of 

Eqs.(13) and (14) the effective pure state can be written generally and con- 
veniently as 

p,ff = 2^'^Do = \E + (2"-i - 1) X T-Huhz-Inz 

n 

- E {exp(-zf7fc,.)H^I^W^exp(2f7,,,) 
fe=i 

x(2"-i7i,72,...7„,)exp(-zf7fe^)W^DW^exp(if7ik,)} 

- E {exp(-if X 2HkxIixImx)WDWeM-iUhx + hx + Imx)\ 

x{2^-^hzhz...Inz) exp[if (7fe, + Iix + Imx)\WDWe-KY>{il X 2^1kxllxlmx)} 

E {exp(-if X 2'lkxIlxImxW<ix)WDW 

q>p>m>l>k=l 

X exp[if (7fe^ + Ii^ + Imx + Ipx + ^<zx)](2""^7i;, 

l2z---Inz) 

X exp[-if (7fc^ + 7,^ + 7^^ + 7p^ + 7<,^)]W^Diyexp(if x 2^1kxlixlmxlpxlqx)} 



It follows from the above expression that one may need to perform (2"^^ + 1) 
different experiments by starting from the density operator, i.e., the longi- 
tudinal n-spin order component a{0) = 2^^^^I\zl2z---Inz in order to prepare 
experimentally the effective pure state 2"~^7?o — \E with n qubits. This 
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shows that the experment number to prepare the n-qubit effective pure state 
by the present method is about half number [2^ — 1) required by the exhaus- 
tive averaging method in Ref.[23]. 

3. Discussion 

It has been shown that the effective pure state can be prepared conve- 
niently on the basis of Eq.(13) when the starting reduced density operator 
is chosen suitably. Experimentally, the effective pure state is prepared gen- 
erally from the starting density operator, that is, the longitudinal n-spin 
order component, a{0) — 2^~^Iizl2z---Inz in a unified and systematical form, 
while the latter can be obtained conveniently from the equilibrium magne- 
tization in a coupled n-spin (1=1/2) system. Actually, the starting density 
operator (t(0) = 2'^~^Iizl2z---Inz can be obtained easily from the equilibrium 
magnetization aueq oc Ikz of any kih. spin in the coupled spin (1=1/2) sys- 
tem by performing a sequence of one-qubit quantum gate operations and 
the two-qubit diagonal gate operations Gki{Xki) — exp{—iXki2IkzIiz) [26], for 
example, 

(^leq liz ^■llx-i2z ^-'lz-'2z ^^lz^2z^3z 

T T TT r 

ori-2r T T T^J-n-ly-l-nz 2 "-^i/ on-1 t T T —r^(f\\ 
^ J-lzJ-2z---J-n-lz ^ J-lzJ-2z---J-nz — O^^UJ, 

where the unitary transformation exp(— i^^P)A exp(?^P) = 5 is denoted 

ep 

briefly as A — B. In the high-temperature approximation the thermal equi- 
librium state of an n-spin system can be written as a^q oc Yl,k ^kz- Then 
(Tkeq oc Ikz Can be obtained by first applying a selective 90° pulse only to 
the spin k, then a hard QO!.^. pulse to all the spins of the system, and then 
applying a magnetic field gradient to the system to cancel all the trans- 
verse magnetizations Iiy (/ 7^ k). The equilibrium magnetization of any 
one of the spins in the spin system, e.g., the kih spin, usually can be en- 
hanced easily by various polarization transfer processes [11, 24, 25]. Then 
the present method is a simple scheme to make use of various polarization 
transfer techniques to scale up the qubit size on an NMR quantum computer 
[11, 21]. Because the equilibrium magnetization is inverse proportional to 
the exponential factor 2" when the high-temperature approximation holds 
in the coupled n-spin system it is easy to find that the effective pure state 
prepared by the method has a signal-to-noise ratio inverse proportional to 
the factor 2" in each experiment on average, and it is well-known that this 
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is a general disadvantage of ensemble quantum computation [11, 12, 23]. 

The present preparation of the effective pure state involves the unitary op- 
erations including the n-qubit Walsh-Hadamard transform W, the n-qubit 
diffusion transform D [5], and the m-body (1 < m < n) elementary prop- 
agators exp{±i92"^~^ Ikixh2x- ■ -hmx) built up with the m-body interaction 
2"^-'4ix42x-4^x (1 < A;^ < n, and i = l,2,...,m) [26, 27]. Each of all 
these unitary operations is readily decomposed completely into a sequence 
of one-qubit quantum gates and the two-qubit diagonal quantum gates [26]. 
Therefore, the present method to prepare the effective pure states could be 
programmed and performed conveniently on an NMR quantum computer. 

Acknowledgement 

This work was supported by the NSFC general project with grant number 
19974064. 

References 

1. R.Feynman, Int.J.Theoret.Phys. 21, 467 (1982); Found. Phys. 16, 507 

(1986) 

2. D.Deutsch, Proc.Roy.Soc.Lond. A 400, 97 (1985) 

3. P.W.Shor, Proceedings of the 35th Annual Symposium on Foundations of 
Computer Science, edited by S.Goldwasser (IEEE Computer Society Press, 
Los Alamitos, CA, 1994), p. 124; SIAM J.Comput. 26, 1484 (1997) 

4. D.Deutsch and R.Jozsa, Proc.Roy.Soc.Lond. A 439, 553 (1992) 

5. L.K.Grover, Phys.Rev.Lett. 79, 325 (1997) 

6. J.I.Cirac and P.ZoUer, Phys.Rev.Lett. 74, 4091 (1995) 

7. Q.A.Turchette, C.Hood, W.Lange, H.Mabushi, and H.J.Kimble, 
Phys.Rev.Lett. 75, 4710 (1995) 

8. D.P.DiVincenzo, Science 270, 255 (1995) 

9. Yu.Makhlin, G.Schon, and A.Shnirman, Nature 398, 305 (1999) 

10. S.Lloyd, Science 261, 1569 (1993) 

11. N.Gershenfeld and I.L.Chuang, Science 275, 350 (1997) 

12. D.G.Cory, A.F.Fahmy, and T.F.Havel, Proc.Natl.Acad.Sci.USA 94, 1634 

(1997) 

13. I.L.Chuang, L.M.K.Vandersypen, X.Zhou, D.W.Leung, and S.Lloyd, Na- 
ture 393, 143 (1998) 

14. I.L.Chuang, N.A.Gershenfeld, and M.Kubinec, Phys.Rev.Lett. 80, 3408 

(1998) 

15. J.A.Jones, M.Mosca, and R.H.Hansen, Nature 393, 344 (1998) 



11 



16. D.G.Cory, M.D.Price, W.Maas, E.Knill, R.Laflamme, W.H.Zurek, 
T.F.Havel, and S.S.Somaroo, Phys.Rev.Lett. 81, 2152 (1998) 

17. Z.L.Madi, R.Bruschweiler, and R.R.Ernst, J.Chem.Phys. 109, 10603 
(1998) 

18. N.Linden, H.Barjat, and R.Freeman, Chem.Phys.Lett. 296, 61 (1998) 

19. M.A.Nielsen, E.Knill, and R.Laflamme, Nature 396, 52 (1998) 

20. S.S.Somaroo, C.H.Tseng, T.F.Havel, R.Laflamme, and D.G.Cory, 
Phys.Rev.Lett. 82, 5381 (1999) 

21. I.L.Chuang, N.Gershenfeld, M.G.Kubinec, and D.W.Leung, Proc.Roy.Soc. 
Lond. A 454, 447 (1998) 

22. D.G.Cory, M.D.Price, and T.F.Havel, Physica D 120, 82 (1998) 

23. E.Knill, I.L.Chuang, and R.Laflamme, Phys.Rev. A 57, 3348 (1998) 

24. R.Freeman, A Handbook of Nuclear Magnetic Resonance (Longman, Har- 
law, 1987) 

25. R.R.Ernst, G.Bodenhausen, and A.Wokaun, Principles of Nuclear Mag- 
netic Resonance in One and Two Dimensions (Oxford University Press, Ox- 
ford, 1987) 

26. X.Miao, |http://xxx.lanl.gov/abs/quant-ph/0003068| (2000) 

27. X.Miao, Molec.Phys. 98, 625 (2000) 



12 



